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1 Introduction 

Quantum Entanglement (QE)[T][2l[3l|4], the most surprising non-classical property of a 
quantum system, plays a key role in quantum information and quantum computation 
processing. Similarly, topological entanglement (TE) [5] is described in terms of link di- 
agrams and via the Artin braid group. There are natural relationships between QE 
and TE[6][7]- Kauffman and his co-workers have explored the role of the unitary solu- 
tions to the Yang-Baxter EquationfYBE) [8ll9lll0jin quantum computation. According 
to their theories, the unitary Yang-Baxter R matrices are both universal for quantum 
computation and are also solutions to the condition for topological braiding. This mo- 
tivates a novel way to study YBE(as well as braid relations) 11 . 12 , 13. 14. 15 .161117] . A 
set of size 4x4 universal quantum gates are constructed in terms of unitary R ma- 
trices, for example, the CNOT gate[6], DCNOT gate(i.e.. Double CNOT gate)[Tg. By 
means of universal R matrix, entanglement swapping and Yang-Baxter Hamiltonian 
are investigated in Ref. [TS]. In Ref.[15], Chen et al. point out that all pure two-qudit 
entangled states can be achieved via a universal Yang-Baxter R matrix assisted by 
local unitary transformation. Later on, the geometric properties of this Yang-Baxter 
system is studied in Ref.[19]. 

Temperley-Lieb algebras(TLA) grew out of a study of solvable lattice models in 
two-dimensional Statistical Mechanics |20j and is related to link and knot invariants [21]. 
a recent study [22^ shows that TLA is found to present a suitable mathematical frame- 
work for describing quantum teleportation, entanglement swapping, universal quantum 
computation and quantum computation flow. Additionally, the systems of qutrits or 
more generally qudits are more powerful than the systems of qubits habitually used in 
quantum computer 23, 24, 25, 26, 27, 28, 29^. Due to the importance of TLA in quantum 
information processing, we find matrix realizations of TLA in high dimension in this 
paper. Consequently, by means of Yang-Baxterization approach, a family of universal 

X v? R matrices associated with TLA can be constructed. 

This paper is organized as follows: In SecO we recall the method of construct- 
ing matrix realizations of TLA which is given by P.P.Kulish. Then we present a 
method of constructing x v? matrix realizations of TLA with n"^ nonzero matrix 
elements. In SecO a unitary v? x Yang-Baxter R matrix is constructed via Yang- 
Baxterization[30| acting on the x n? matrix realizations of TLA. In Secl4l when 
n=3, we investigate the entanglement properties of _R(^, gi, g2)-matrix. We show that 
arbitrary degree of entanglement for two-qutrit entangled states can be generated via 
the unitary R{6,qi,q2)-raa.triK acting on the standard basis. Then we can construct a 
Hamiltonian from the unitary R{6, gi, g2)-niatrix. Furthermore, the Berry phase of the 
system is investigated, and the results show that the Berry phase of this system can be 
interpreted under the framework of SU(2) algebra. This result is consistent with that 
given in Ref. [19]. 



2 An extended method of constructing realizations of TLA 

In this paper,the matrix realization of TLA U-matrix and YBE solution _R— matrix are 
v? X matrices acting on the tensor product space V x V, where V is a n— dimensional 
vector space. As U and R act on the tensor product Vi x Vi+i, we denote them by Ui 
and Ri, respectively. 
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We first briefly review the theory of TLA [20]. For each natural number m, the TLA 
TLm{d) is generated by {7, Ui,U2 ■ ■ ■ Um~l} with the TLA relations: 

Uf =dUi 1 < i < m - 1 

UiUi±iU^^U^l<i<m (1) 

U,Uj=UjU, \i-j\>2 

where the notation Ui = Ui^i^i is used. The Ui represents li (8) I2 ® I3 ■ ■ ■ ® Ij-i ® f^® 
li+2 <8> ■ • ■ iX) Im , and Ij represents the unit matrix in the jth space Vj . In topology, 
the parameter d corresponds to a single loop "O"- In addition, the TLA is easily 
understood in terms of knot diagrams in Ref.[6]. 

In Ref. [31], P.P.Kulish et al. showed a method of constructing matrix realizations 
of TLA. Let us review it briefly. For a given invertible nxn matrix A, a n^ x n^ matrix 
solution can be constructed in terms of A and A"^ with [/"J" = A^(y4~^)^. Hereafter, 
(7°^ denotes C/^b cd ^-nd Af^ denotes A^ f, with a,b,c,d = 0, 1, 2, • • • , n — 1. One can 
verify that U is a matrix realization of TLA. Let Tr{M) denote the trace of matrix M, 
and M'^ denote the transpose of matrix M. In terms of A and A~^, the single loop d 
can be determined by d = Tr{A^ A~^). By means of this method, many realizations 
of TLA can be constructed. For example, we set 

Then a 4 x 4 matrix realization of TLA can be constructed as follows 

/ 1 g 
I 

\q^^ 1 

However, Not all solutions can be constructed with this method (for example, the 
solution associated with eight vertex model can't be constructed). In the following, we 
will introduce an extended method of constructing matrix realizations of TLA. 

In order to construct a matrix realization of TLA, we introduce two nxn invertible 
matrices A and B. We assume that matrix U can be constructed as t/^J' = A^B^- 
Substituting this relation into TLA relations ([T}, the limited conditions for A and B 
can be derived. The relation C/^ = dU yields d = Tr{B^ A). Then U is a realization of 
TLA relations ((T} if and only if A and B respect the following conditions 

(BAfiAB) = {AB){BAf = J„xn. (3) 

Where Inxn represents the unit matrix in n dimension. In particular, if we take 
B = A~^, we note that the condition ([3]) is obviously satisfled. Thus we re-obtain 
P.P.Kulish's method of constructing matrix realizations of TLA. 

In order to find matrices A and B satisfying these conditions, a special matrix 
structure is adopted, that is, each row and each column has one matrix element, and 
the matrix element locations are on the main diagonal symmetry. In addition, A and 
B satisfy the relation B? — (A^)^^ for the non-vanishing entries. In addition, we can 
verify that relations A B = B'^A = AB'^ = BA^ ^ Inxn hold. In this case, Eq.© 
hold. Then we obtain a matrix realization of TLA. One can verify that Tr{I 
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In fact, we can select n matrices which satisfy these conditions, and all their matrix 
elements occupy different locations. Let i denote the ith matrix. Namely, the non- 
vanishing matrix elements of are {A'^^)}-2, (^^'^)f-3, ■•■ , (^^'^)o"^ 
{A^^^)\_^i, •■■ , {A^^^)^~^ . For example, if n=4 and i—2, the non-vanishing matrix 
elements of A^^) are [A^^^)^, {A'^^'')},, {A^^^)l and {A^^^)^. There are n invertible 
matrix B^^^ which is determined by [i?^*-']^ = for non-varnishing matrix 
elements. In terms of matrices and B^^\ a matrix representation of TLA can be 
constructed as [(7^'^]°^ = [B^*)]^. By means of these n matrix realizations of 
TLA U^^\ we can construct a x matrix realization of TLA with nonzero 
matrix elements. Taking the summation of these n matrices, we can obtain a combined 
matrix 

n 

u^±,j:ui-K (4) 

i—1 

In fact, the realizations of TLA matrix of this form can be represented in terms of 
Dirac's "bra" and "ket". This notation will appear in Sec[2] If we substitute Eq.(|4]) 
into Eqs.lll, the first equation in Eqs.((T} is satisfied automatically, and one can check 
that d = v^n. Then the other two relations are satisfied by the following limiting 
conditions. 



n 

^(^0)bW)(bO)^W)T^O^^^, (5) 

Where i k and i,k — 1, 2, ■ • • n, and 0„ denotes nxn matrix with all matrix elements 
are zero. This limiting condition together with the special matrix structure are used 
to determine U matrix. Two examples are shown to illustrate the application of this 
method in detail. 



2.1 Example I: The case n = 2 



The simplest example which illustrates the method is the case n—2. According to the 
above analysis, when n = 2, we choose two sets of 2 x 2 invertible matrices as follows, 

fid) ^ 



ai 
6i 



b 



(6) 



02 
&2 



a. 







Where and 6^ are parameters which will be determined by the conditions in Eq.©. 
Then two U matrices can be obtained as follows (we choose {jOO), |01), \W), |11)} as 
standard basis), 

o\ 

02^2"^ 

1 
yo 0/ 



[/(I) 



1 OOaibj;^^ 



1 



\a-^bi 



t/(2) 







a, ^62 



(7) 
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The trace of these two solutions is 2 (i.e., di = d2 = 2). 

In order to obtain a solution associated with eight-vertex model, we consider the 
combinatorial structure of 

[/(I) 

and 

{/(2) 

. The combinatorial form reads 

v2 

If we substitute this relation into Eqs.([5]). Then we can derive a strong limiting con- 
dition 026^^ = ei {e — ±). Let M* denote complex conjugation of matrix M . We can 
introduce a new parameter q with q — aib^^, which is complex and has norm l{i.e. 
q* — q^^)- Then a eight- vertex matrix representation with d = \/2 is obtained as 
follows. 



(8) 



Let 

|^i) = -^(100) + g-llll)), 

|V2) = ^(|01)-e*|10)). 
Then, in terms of "bra" and "ket", the U matrix takes the following form 

(7 = %/2(|Vl){Vl| + |^2){^2|). 

This realization of TLA is associated with eight vertex model '9 . And this solution has 
been applied to many fields, such as topological quantum computation[32) and two 
dimensional representation of YBE[16). 



2.2 Example II: The case n = 3 

Let and B^*) (i=l,2,3) are three sets of 3 x 3 matrices with standard basis(j.e. 
|0),|1),|2)). We set 




^(1)^ 



^ 5(2) = (9) 



^(3) 



Where Oj, 6j and are also undetermined parameters. Thus we note that the relation 
Eq.Qis clearly satisfied. If we choose {|00), |01), |02), |10), |12), |20), |21), |22)} as 
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standard basis, then we can obtain three sets of 3^ x 3^ matrices and 
U^^\ In this case, their single loop di = 3(i=l, 2, 3). Then the combined form of U 
matrix U = (f/^^^ + U'-'^^ + U'^^^)/V3. Substituting this combined form into Eqs.©, 
the undetermined parameters follows from the limited conditions, 

(10) 

«lCi ^ = 1, a2C^ ^ = ujq2, a^c^ ^ = qi. 

Where qi = e^"^^ and u satisfies the relation lo^ + lo + 1 — {i.e., cu — e'"'~). On the 
standard basis U has the matrix form 



1 














uqi 





<7l 


\ 





1 





Ul 














Loq2 








1 





91 

92 







1 











1_ 





1 











92 








92 

91 





1 





92 
91 








1 














1 





J_ 

u 











1 





91 

92 





1 








1 

91 
















LU 





1 





1 

W(J2 





1 

92 














1 / 



The single loop of this solution is d = %/3- In fact, we can introduce three sets maximally 
entangled states as 

|Vi) = ^(|02) + gig2-i|ll) + |20)) 
IV2) = -^(|01) +a;-lilO) +u;-V'|22)) 

m = -^m+Lo-\^'\i2) + gr'i2i)). 

In terms of these maximally entangled states, the U matrix can be written in a 
elegant form 

U = %/3(|V'l){Vl| + IV'2){V'2| + |V'3){V'3l). 



2.3 Remarks 



We close this section with some remarks. Juan Ospina made a Mathematica imple- 
mentation of this method, and the results in this article were re- obtained [33]. When 
n = 2, the solution ([8]) has been discussed in many works. As we all know, when n = 3, 
the solution (|lip is not discussed. We note that the solutions([8j and are Hermitian 
matrices (i.e., L*'^ = [/)(This fact will be used in the process of Yang-Baxterization 
approach) . 
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3 Yang-Baxterization of U matrix 

In order to discuss tiie non- maximally entangled states, the author In Ref. [13], the 
unitary R matrix has been introduced in Ref. [13]. To make the paper self-contained, 
we briefly review it in the following. In this work, we utilize the so caUed relativistic 
Yang-Baxter Equation(YBE)[l6]. The relativistic YBE reads, 

^ ^ / U V \ ^ ^ ^ f u ~\~ v \ ^ 

Rr{u)Rr+l T— 5^ Rr{v) = Rr+l{v)R; 7—5^ Rr+li^) (12) 

where Ri represents li ® I2 ® I3 (8) • • • ® li_i (8) -R ® ® ■ ■ ■ ^ Im. The variables u 
and i; are caUed as the spectral parameters. The /3 is a constant with /3~^ = ic (c is 
light velocity). 

Let the unitary R{u) matrix take the form 

R{u)^F{u)[Inxn+G{u)U]. (13) 

Where the functions F{u) and G{u) are to be determined. Substituting Ea. (|13|l into 
Ea. (|12p . we obtain the relation 

G{u) + G{v) + G ( .""to ] [G{u)G{v) -1] + V^G{u)G{v) = (14) 
\ 1 -|- p^'uv j 

Following Hu et aZ. [16], we set 

G(u 



aj3u 



bp'^u^ + c/3u -I- d' 

If we substitute it to the relation Ea. (ll4|l . Then we obtain equations for undetermined 

parameters a, b, c and d, 

' c? + Vnac + -f 3fed -I- = 
< ^/na -I- 2c = 
, 6 = d 

After some algebra, a solution of G(u) is obtained as follows 

G(n) = . (15) 

y/4~^{(3'^u'2 - 2^yn/{A-n)iel3u + 1) 

We note that n 7^ 4. The case d — y/2 has been discussed in Ref. [16]. In addition, the 
unitary relation R} {u)R(u) = R{u)R^ (u) = Inxn leads to the relation F* {u)F{u) — 1 
and G{u) + G* {u) + ^/nG{u)G* (u) = 0, where * denotes complex conjugation. Consider 
these relations, one can introduce a new variable 9 with G{u) — (e^^*^ — l)/^/n, which 
equivalent to the relation 



/3^u^ + 2^n/(4 - n)ie(3u + 1 
^92^2 - 2^yn/{4-n)iel3u + 1 



■ e 



-2ie 
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We set F{u) — e'^ with 9 is real. In terms of the new variable, we rewrite the Yang- 
Baxter matrix in a new form 

Ufa \ iSj 2isine 

R(0,qi,q2)^e Inxn l^U. (16) 

The case of n—2 has been discussed in Ref.(|16j). If n—3, on the standard basis the 
unitary solution of R matrix is 

















^gqi 





991 


\ 





/ 





tog 














Logq2 








/ 





9^ 





9 











g. 





/ 



























/ 





y 91 








a 














/ 





a. 











9 





g- 

y 92 





/ 








9 
91 














ujg 





/ 





V 


a 

ujq2 





a 

12 














/ / 



Where / = f{e) = (e"^^ + 2e^^)/V3 and g = g{e) = (e^*^ - e*^)/\/3. 



4 Entanglement and Hamiltonian 

By Brylinskis theorem [34] . a 4 x 4 Yang-Baxter R matrix is universal for quantum 
computation, if and only if this Yang-Baxter R matrix can generate entangled states 
from separable states. The proof of universality for x Yang-Baxter matrix is 
presented in Ref.[15). Via a unitary universal Yang-Baxter R matrix acting on the 
standard basis, one can obtain a set of entangled states. For example, if one lets R{9) 
act on the separable state \lm)[i.e., \l) ® |m)), this yields the following family of states 
\i^)lm ~ y^fj— n" ^ R\iyi\l'fn){KfTi ~ 0, 1, ■ ■ ■ ,n — 1). These unitary matrices may be 
universal for quantum computation, hence they can entangle states. The case n=2 has 
been discussed in Ref.[13). 

Hereafter we focus on the case n=3. For example, if 1=0 and m=Q, then i^)oo = 
(/jOO) + u)~^ gq^^\12) -|- (;gj~^|21))/3. By means of negativity, we study these entangled 
states. It should be noted that the negativity criterion is necessary and sufficient only 
for 2(g) 2 and 2(8)3 quantum systems. However, negativity is well-defined for calculation, 
and it has been widely applied to evaluation of entanglement [35II36II37] . The negativity 
criterion for two qutrits is given by 

AA(p) = ^^^, (18) 

where \\p^^ ||i denotes the trace norm of p^^ , p^^ denotes the partial transpose of the 
bipartite state p. The N{p) corresponds to the absolute value of the sum of negative 
eigenvalues of p^^ , and negativity vanishes for unentangled states. Then negativity of 
the state |'(/')oo yields 

N{e) = -{sin^e +\shie\yJl + Scos'^e). (19) 
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If I.9I = K = e' 3 ), then the state |V')oo becomes the maximally entangled state 

for two qutrits 

Wno = -^(e^^lOO) -ia;-lg-lll2) - iq^^\21)). 

In general, the unitary Yang-Baxter matrix R{9) acts on the basis {|00), |01), |02), 
|10>, |11), |12), |20), |21), |22)}, we obtain the same range of negativity as EqUSl. It is 
easy to check that the negativity ranges from to 1 when the parameter runs from 
to TT. But for £ [0,7r], the negativity is not a monotonic function of 9. And when 
9 — 7r/3, R{9) generate nine complete and orthogonal maximally entangled states for 
two qutrits. 

In fact, we can introduce a unitary transformation Y = Yi igi ¥2- Yi and Y2 take 
the form 

/ c*^ \ / e-*^ \ 

Fi = 1 Y2 = 1 1 . 
\ e"*Tf J \ e"*Tf / 

By means of this local transformation, the universal R matrix (|17|l is locally equivalent 
to R matrix in Ref.[19). 

A Hamiltonian of the Yang-Baxter system can be constructed from the R{9, ipi , 'P2)- 
matrix. As shown in Ref.[16j. the Hamiltonian is obtained through the Schrodinger 
evolution of the entangled states. Let the parameters i^i be time- dependent as tpi — ujit. 
The Hamiltonian is 

H = ^n^-^^^^^^R\9,^,,^2). (20) 

This Hamiltonian is equivalent to the Hamiltonian in Ref. [19], so one can obtain the 
same results as in Ref. 19 . The Berry phase of this system can also be explained in the 
framework of SU(2) algebra. The Berry phase can be explained as solid angle which is 
expanded in the parameter space. We will not discuss this in detail in this paper. But 
we should note that the meaning of the parameter 9 is different. The 9 in Ref.f|19j) 
arises from trigonometrical parameterization, and the 9 in this work arises from the 
relativistic rational parameter. 

5 Summary 

In this paper, we present a method of constructing n? x matrix realization of 
TLA. This matrix realization of TLA has n"^ nonzero matrix elements. Applying Yang- 
Baxterization approach to the matrix realization of TLA, one can obtain & x v? 
Yang-Baxter R matrix. When a Yang-Baxter R matrix acts on the standard basis, one 
can obtain a family of entangled states. Yang-Baxter R matrix is universal for quantum 
computation. 

We believe that this family of Yang-Baxter R matrices associated with U matrices 
will be applied in quantum information, quantum computation and so on. We will 
investigate these applications in subsequent papers. 

Acknowledgements The authors gratefully acknowledge Juan Ospina for helpful comments 
on this paper. Special thanks to the first referee for his advice and criticism on our manuscript. 
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